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Abstract: We argue that the whole or the part of the phase diagrams of QCD and QCD- 
hke theories should be universal in the large- A'^c limit through the orbifold equivalence. The 
whole phase diagrams, including the chiral phase transitions and the BEC-BCS crossover 
regions, are identical between SU(A'^c) QCD at finite isospin chemical potential and S0(2A'^c) 
and Sp(2A'^c) gauge theories at finite baryon chemical potential. Outside the BEC-BCS 
crossover regions in these theories, the phase diagrams are also identical to that of SU(A'^c) 
QCD at finite baryon chemical potential. We give examples of the universality in some 
solvable cases: (i) QCD and QCD-like theories at asymptotically high density where the 
controlled weak-coupling calculations are possible, (ii) chiral random matrix theories of 
different universality classes, which are solvable large- (large volume) matrix models of 
QCD. Our results strongly suggest that the chiral phase transition and the QCD critical 
point at finite baryon chemical potential can be studied using sign-free theories, such as 
QCD at finite isospin chemical potential, in lattice simulations. 
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1. Introduction 



Unraveling the properties of quantum chromodynamics (QCD) at finite temperature T 
and finite baryon chemical potential /x_b is essential to understand various phenomena from 
ultrarelativistic heavy-ion collisions, the early Universe, supernova explosion, and neutron 
stars, to possible quark stars (for a recent review, see, e.g., [^). A lot of progress have 
been made theoretically by the first-principles lattice simulations in QCD at high T and 
sufficiently small /U^ in connection with the experimental investigation of the quark-gluon 
plasma in the ultrarelativistic heavy-ion collisions at RHIC and LHC. On the other hand, 
the lattice technique is not available at finite because of the sign problem: the fermion 
determinant is no longer real at finite /i^, and the Monte Carlo approach based on the 
importance sampling does not work. This is why QCD at finite hb has not been fully 
understood, such as the precise location of the QCD critical point(s) (for a review, see |2|) 
and the realization of the color superconductivity at intermediate ^b (for a review, see ||3|) 
relevant to the physics of neutron stars, etc. 

Still there are a class of QCD-like theories which are free from the sign problem. These 
theories intensively studied so far include two-color QCD with even numbers of fundamental 
flavors [|, I, any-color SU(7Vc) QCD with adjoint fermions [|], and SU(7Vc) QCD at finite 
isospin chemical potential /U/ (^,0]. In addition to the chiral phase transition, these theories 
exhibit the Bardeen-Cooper-Schrieffer (BCS) pairing of quarks at large chemical potential 
due to the same mechanism as the color superconductivity in real QCD However, it is 
not clear, at the quantitative level, or even at the qualitative level, how the phase diagrams 
of these theories are related to each other, and more importantly, to that of real QCD at 
finite ^b- There are also other theories which are free from sign problem at finite ^b but 
have not been well studied: S0(2A'c) and Sp(2A^c) gauge theories with Nf fundamental 
Dirac fermions.^ Apparently, there is no a priori reason why these theories capture the 
physics of SU(A'^c) QCD at finite chemical potential either, since the structures of gauge 
groups are different. 

In this paper, we argue that the whole or the part of the phase diagrams of QCD and 
QCD-like theories should be universal in the large- A'^c limit via the orbifold equivalence.^ 
The relations between QCD and QCD-like theories are summarized in Fig. |^. The whole 
phase diagrams described by the chiral condensate and the superfluid gap should be quan- 
titatively identical between SU(A^c) QCD at finite /i/ and S0(2A^c) and Sp(2A^c) gauge 
theories at finite ^b-^ As common properties of these theories, the Bose-Einstein conden- 

^The symplectic group is defined as Sp(2Ai'c) = {g £ U(2A^c)|<7^ Jeff = Jc}, where Jc is the antisymmetric 
matrix defined in ( |2.10[ ). The dimension of Sp(2iVc) is Nc{2Nc + 1). 

^The universahty of the Wilson loops between SU(A^c), SO(2A'^c), and Sp(2A'^c) in pure gauge theories in 
the large- A^c limit was pointed out long time ago in |^] , even before the notion of the orbifold equivalence was 
found. Probably the first paper which considered an equivalence of the phase diagrams of QCD-like theories 
is 1^, in which SU(A'^c) gauge theory with two index fermion representations [adjoint, (anti)symmetric, and 
bifundamental] on the compact space x S''^ as a function of volume was studied. Our paper is the first, 
to our knowledge, to argue the large- A'^c universality of the phase diagrams of QCD and sign-free S0(2A'^c) 
and Sp(2A'^c) gauge theories with fundamental fermions at finite chemical potential on physically relevant 
space R*. 

^We note that the deconfinement phase transition characterized by the Polyakov loop should also be 
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sate (BEC) of the Nambu-Goldstone (NG) modes at small chemical potential, as well as 
the BCS-type pairing at large chemical potential, appear in the phase diagrams. Since the 
BEC and BCS pairings carry the same quantum numbers and break the same symmetry, 
the BEC and BCS regions should be continuously connected without any phase transition. 
This is the BEC-BCS crossover similar to that in nonrelativistic condensed-matter systems 
pO| , ^ and the continuity between the hadron phase and the color superconducting 
phase (hadron-quark continuity) in three-flavor QCD at flnite fiB [0, 0, The univer- 
sality of the phase diagrams means that both the chiral phase transitions and the BEC-BCS 
crossover regions should appear at the same coordinates in the T-fi plane independently of 
the theories in spite of the different symmetry breaking patterns. 

We also argue that the phase diagrams of these theories outside the BEC-BCS crossover 
regions should be identical to that of SU(A'c) QCD at finite ^b, which is most relevant in 
reality. In particular, the magnitude and temperature dependence of the chiral condensate 
in QCD at finite hb should be exactly the same as those in QCD at flnite /i/. Since the 
latter can be obtained by dropping the complex phase of the fermion determinant in the 
former for even number of flavors, this suggests that the phase-quenched approximation 
for the chiral condensate is exact in this region. Actually, this phenomenon has already 



been observed in the lattice QCD simulations |19| and the model calculations, such as the 



chiral random matrix model pO], the Nambu~Jona-Lasinio model |21, p^], and the hadron 



resonance gas model [23|, though the reason has been unclear so far. We note that our 
arguments based on the orbifold equivalence are model-independent (though it is exact in 
the large- A'^c limit), and thus, would provide a solid theoretical basis of the phase-quenched 
approximation . 



The idea of the orbifold equivalence flrst originates from the string theory |24, 26 1 



and then has been generalized to the quantum fleld theory without reference to the string 



theory |27, 29, |3^, ^ 33, 34]. The procedure of orbifolding (or orbifold projection) 
is as follows. Identify a discrete global symmetry of the original theory (the "parent" 
theory) . Eliminate all the degrees of freedom in the parent theory which are not invariant 
under the discrete symmetry. This gives a new theory (the "daughter" theory). Then a 
class of correlation functions and observables are shown to be identical between the parent 
and daughter theories in the large- Ac limit, as long as the symmetry used to make the 



projection is not broken [32|. For example, as recently proposed |3^, SU(Ac) QCD 



with Nf quarks at flnite fiB can be obtained from the S0(2Ac) gauge theory with Nj 



identical between these theories. Because the quark chemical potential fi does not affect the gauge dynamics 
as long as = 0{Nji!) in the large- A^c counting, the deconfinement temperature is independent of fi at the 
leading order; the universality at nonzero /i follows once that at ^ = is provided M. For this reason, 



except Sec. 3.2 where we discuss the 1/Nc correction to the deconfinement temperature, we shall mostly 
concentrate on the phases related to the chiral and flavor dynamics which presumably depend on at the 
leading order of 1/Nc expansion. 

''The BEC regions in these theories appear in a model-independent way (see the discussion in Sec. 
This is in contrast to QCD at finite fj.B where the BEC of diquark pairing can appear depending on the 
details of the models Jl^, p^ . Throughout this paper, we assume the crossover between the BEC and BCS 
regions similarly to [[^, [f8[ though our result of the universality of the phase diagrams does not rely on this 
assumption. 
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S0(2A/c) or Sp(2A/c) YM with jLve 



Equivalent outside 
the BEC-BCS 
crossover region 



Equivalent 
everywhere 



SU(A/c) QCD with [jb 



SU(A/c) QCD with [ji 



Equivalent outside the 
BEC-BCS crossover region 

Figure 1: Relations between SU(iVc) QCD at finite (/i/) and S0(2A^c) and Sp(2iVc) Yang-Mills 
(YM) theories at finite ^b. SU(7Vc) QCD at finite fj.i can be obtained from S0(2iVc) or Sp(27Vc) 
gauge theory at finite fiB through the orbifold projection in the whole phase diagram, while QCD 
at finite /is can be obtained outside the BEC-BCS crossover region of these theories. As a result, 
QCD at finite //_b is equivalent to QCD at finite fij outside the BEC-BCS crossover region of the 
latter. 



fundamental fermions at finite hb through an orbifold projection. Here it is this condition, 
"outside the BEC-BCS region" mentioned above, that U(l) baryon number symmetry of 
the parent S0(2A'^c) and Sp(2A''c) gauge theories, which is used for the projection, is not 
broken spontaneously. 

To be precise, the orbifold equivalence has been proven to all orders in the perturba- 
tion theory but not nonperturbatively, except certain QCD-like theories containing adjoint 
scalar or fermion matter at //s = /i/ = [31| and supersymmetric QCD at finite fiB or m 
in a holographic setup |^]. In this paper, we provide new evidence for the nonperturbative 
orbifold equivalence within real QCD and QCD-like theories with fundamental fermions. 
At sufficiently large chemical potential, nonperturbative observables, such as the super- 
fluid gap and the diquark condensate, can be computed using the controlled weak-coupling 
calculations owing to the asymptotic freedom. Then we can explicitly compute the l/N^ 
corrections and demonstrate that the equivalence is not only exact in the large-A'^c limit 
but is well satisfied even for Nc = 3. We also apply the idea of the orbifold equivalence to 
the chiral random matrix theory, a solvable large- (large volume) matrix model of QCD 
(for reviews, see |]38| , ^]). We verify the nonperturbative orbifold equivalence between the 
random matrix theories of different universality classes by computing the effective poten- 
tials in the N ^ oo limit (thermodynamic limit). Our calculations and arguments, though 
do not constitute a complete proof, provide overwhelming evidence for the nonperturbative 
equivalence in real QCD at finite chemical potential. 

The paper is organized as follows. In Sec. |2|, we study the properties of SU(Ac), 
S0(2A^c), and Sp(2Ac) gauge theories at finite chemical potential. In Sec. ^, after a brief 
review of the orbifold equivalence in the large- Ac QCD, we construct the orbifold projec- 
tions between QCD and QCD-like theories. We then argue that the phase diagrams of these 
theories should be universal. Also we compute the 1/Ac corrections at asymptotically high 
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density. In Sec. we construct the orbifold projections between the chiral random matrix 
theories of the different universahty classes. We then expUcitly check the nonperturbative 
orbifold equivalence. Section ^ is devoted to conclusion and discussion. 

2. Classification of QCD and QCD-Iike theories 

In sections ^ and |3|, we consider the following classes of theories: SU(A''c > 3), S0(2A'^c); 
and Sp(2A^c) Yang-Mills theories with Nf fundamental Dirac fermions. For SU(A'^c) gauge 
theory, we consider finite baryon chemical potential /U^ or finite isospin chemical potential 
/X/. For S0(2A''c) and Sp(2A''c) gauge theories, we consider finite baryon chemical potential 
H B ■ We can classify these theories by the Dyson index /3 of the Dirac operator which reflects 
the number of independent degrees of freedom per matrix element in corresponding chiral 
random matrix theories, as we will briefly review in Sec. ^: SU(A''c), S0(2A''c), and Sp(2A''c) 
gauge theories correspond to /3 = 2, /3 = 4, and /3 = 1, respectively. 

The Lagrangian of the gauge theories in the Euclidean spacetime is given by 

Nf 

= A t^^/f.)' + E v^/ + HV'f , (2.1) 

where G denotes the gauge group SU(A''c), S0(2A''c), or Sp(2A'^c) and / denotes the flavor 
index. is the field strength of each gauge field = A^^T^ with being the 

generators of each gauge group normalized such that tr{T^T^) = {l/2)6ab- The Dirac 
fermion tp'^ belongs to the fundamental representation of the gauge group G,^ and nif = m 
is the degenerate quark mass.^ The Dirac operator D is defined as 

V = 7^1)^ + m\ (2.2) 
for finite quark chemical potential /u, and 

V = rD,, + \^Jin^T\ (2.3) 

for finite isospin chemical potential /x/ = 2ii when Nj = 2J Here = + iA^ is 
the color covariant derivative. The SU(2) isospin generators Tj are normalized such that 
tr(rjrj) = 25ij. In both cases above, the Dirac operator preserves the chiral symmetry 
{75, = 0. 

In this section, we will investigate the properties (the phase diagram in particular) of 
each gauge theory. For completeness, we repeat some of arguments on SU(A^c) QCD at 
finite /x^ and fij in the literature. The arguments in this section are independent of 
Nc, except that we consider Nc > 3 for SU(A'^c) QCD. 

^For G = SO(2A^c), it is the 2Afc-component vector representation. 

®The degenerate quark mass is not essential in our argument. In order to avoid the sign problem, 
pairwise nif is necessary in SU(A'^c) QCD at finite /ij and Sp{2Nc) gauge theory at finite /is, but is not in 
SO{2Nc) gauge theory at finite /is- See the subsections below. 

'^In this paper, we consider the isospin chemical potential fij for Nf — 2 unless otherwise stated explicitly. 
Our argument can be generalized to any even Nf if we define fij by regarding Nf /2 quarks as "up" with 
the quark chemical potential fi and Nf/2 quarks as "down" with the quark chemical potential — /i. 
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2.1 QCD: SU(iVc > 3) gauge theory at finite ^13 (/? = 2) 

Let us first consider the SU(A''c > 3) gauge theory with fundamental fermions (QCD). In the 
QCD vacuum (// = 0), the Dirac operator V is anti-Hermitian and its eigenvalue i\n defined 
by Vijjn = iXn'4'n is always pure imaginary A^, G M. From the chiral symmetry {75,!?} = 0, 
—i\n is also the eigenvalue when A„ 7^ and the fermion determinant det(P + m) is always 
real and nonnegative. This allows us to use the standard Monte-Carlo simulation based on 
the importance sampling. However, this is no longer true in QCD at finite fiB- Because 
An is generally complex and so is the fermion determinant, the standard Monte-Carlo 
simulation technique fails (the fermion sign problem). In this case, the Dirac operator is 
written in the form of the complex matrix [see ( |4.3D ], and the Dyson index is /3 = 2. 

Although present understanding of the QCD phase diagram at intermediate is 
largely model-dependent the ground state of QCD at low T and at sufficiently large 
/is is expected to be a color superconductor @ based on the controlled weak-coupling 
calculations. At sufficiently large f^B, the physics near the Fermi surface is described by 
the weakly interacting quarks due to the asymptotic freedom. The perturbative one-gluon 
exchange interaction, 

CoGE = -Gii^YTf'^l.f, (2.4) 

is dominant at large fiB with G > 0. The color factor of this interaction in the iptp-chaimel 
is 

where and Tj^ are symmetric and antisymmetric Nc x Nc Hermitian matrices. This is 
attractive in the color antisymmetric channel. According to the Bardeen-Cooper-Schrieffer 
(BCS) mechanism, any infinitesimally small attractive interaction between quarks leads to 
the condensation of quark-quark pairs, the diquark condensate. Because the positive parity 
state is favored by instanton effects the pairing in the spin-parity 0^ channel is 

the most favorable energetically. Considering the Pauli principle, the condensate must be 
flavor antisymmetric. Therefore, the condensate {ip'^ C^^Tata''^) is formed at large /i^, 
where both A and A' denote antisymmetric representations. Since the diquark condensate 
breaks SU(A'^c) color symmetry, this is called the color superconductivity [Q. 

2.2 QCD: SU(iVc > 3) gauge theory at finite {(3 = 2) 

Now we turn to QCD at finite /i/ and /i^ = (see |^] for details). Due to the property 
TiTsTi = — Ts, the Dirac operator ( p^ ) satisfies 

ri75^?75n = (2.6) 

From this relation and the chiral symmetry {75,2?} = 0, if we take an eigenvalue zA„ of 
the Dirac operator (|2.3D , eigenvalues appear in quartet (iA^, — iA„, iA* , — iA* ). Therefore, 
one finds that det[P(/i) + m] > and the standard Monte-Carlo simulation technique can 
be used at finite /ij. 
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When m = /i/ = 0, the Lagrangian has S\J(2)l x SU(2)/j symmetry. If the degenerate 
quark mass m is turned on, the symmetry is exphcitly broken to SV{2)l+r. The isospin 
chemical potential /i/ further breaks SU(2)j;,+/{ down to \J{1) l+r- Note that this symmetry 
is different from \J{1)b- 

Let us consider the zero-temperature ground state at small fij and sufficiently large 
m where the theory is analytically controllable. Unlike the phase structure of QCD at 
finite hb, that of QCD at finite fij (with degenerate quark mass) is rather well understood 
because of the absence of the sign problem and the constraints from the rigorous QCD 
inequalities |42, 4^, 44, ^] [Q, ||, where both properties follow from the relation ( |2.6| ). The 



phase diagram in the T-^u plane is summarized in Fig. For small below the p meson 
mass rup, we can concentrate on the pions at low-energy. When /i/ < m^, no particles 
can be excited so that the ground state is the same as the QCD vacuum at fij = 0. On 
the other hand, when > m-,^, it is favorable to excite vr"*" whose excitation energy is 
"T-TT — ^/ < 0; the ground state turns into the Bose-Einstein condensation (BEC) phase, 
(7r+) / 0, where \J{1)l+r is spontaneously broken down to Z2. 

At sufficiently large /i/, on the other hand, \J{1)l+r symmetry is spontaneously broken 
by the BCS-type diquark pairing. The color factor of the one-gluon exchange interaction 
in the ^-i/'-channel is 

iT!'')a,iT!%s = - ^iT!^UiT!%,, (2.7) 



which is attractive in the color-singlet channel. The spin-parity 0~ channel is favored by 
the instanton effects, and the condensate must be flavor antisymmetric from the Pauli 
principle, (d'y^u) 7^ 0. This is consistent with the requirement from the QCD inequality 
0. Note that, although the mechanism leading to the pairing is similar to the color 
superconductivity in QCD at large hb, this condensate has different quantum numbers 
and does not break SU(A'^c) color symmetry. 

Since the BCS pairing {d^y^u) 7^ at large /i/ has the same quantum numbers and 
breaks the same 11(1)^+^2 symmetry as the BEC (vr"'') 7^ at small it is natural to 
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expect that the BEC and BCS regions are continuously connected without any phase 
transition (the BEC-BCS crossover). At sufficiently high T, the condensate melts away 
and \J{1)l+r symmetry recovers. The critical temperature Tc vanishes at ^ = ?n,^/2 and 
is expected to be an increasing function of fi in both BEC and BCS regions. 

2.3 SO(2iVc) gauge theory at finite = 4) 

We next consider S0(2A'^c) gauge theory at finite fj,B- From the property (A^^)* = —A^^, 
one has the relation: 



C-f5VCj5 = 'D*, or [P,C75i^] = 0, (2.8) 

where C is the charge conjugation matrix satisfying = 1 and Cj^C = — 7^^* with 
all 7 matrices being Hermitian, and K is the complex conjugation operator. The Dirac 



operator can be written in the form of quaternion real matrix [see (4.5)] and the Dyson 
index is /3 = 4, the same universality class as SU(A'^c) gauge theory with adjoint fermions 
Q (see also EH]). From (2.8) and the chiral symmetry {75,2?} = 0, if iXn is one of the 



eigenvalues of V, eigenvalues appear in quartet (iAn, — iAn, iA* , — iA* ). Note that, when 
A„ is real or pure imaginary, this quartet reduces to two sets of doublets {iXn, —iXn) 
with their eigenvectors being linearly independent from the anti-unitary symmetry ( p.^ ) 
||47[ |. Therefore, det[P(^) + m] > and the standard Monte-Carlo simulation technique is 
available at finite /is [B5|, |3^]. 



When m = /is = 0, the Lagrangian (l2J|) has the SU(A^/)l x SU(iV/)R x U(l)s x U(1)a 
symmetry at the classical level at first sight. However, chiral symmetry of the theory is 
enhanced to \J{2Nf) owing to the anti-unitary symmetry (|2.8| ) [p8| , At the quantum 
level, U(1)a C \]{2Nf) is explicitly broken by the axial anomaly and SU(2A'^^) symmetry 
remains. One can actually rewrite the fermionic part of the Lagrangian (|2.lD manifestly 



invariant under S\J{2Nf) using the new variable ^ = (^l, fT2V'/j)"^: 

£f = i^^a^^D^^, (2.9) 

where cj^ = {—i,ak) with the Pauli matrices a^- The chiral symmetry SU(2iVj) is sponta- 
neously broken down to S0(2A'^^) by the formation of the chiral condensate {tpip), leading to 
the 2Nj+Nf-l Nambu-Goldstone (NG) bosons living on the coset space SU(2iV/)/SO(2iV/). 
In contrast to real QCD, there are not only U(l)s neutral NG modes with the quantum 
numbers Ha = Tpj^Paip (just like the usual pions), but also U(1)b charged NG modes with 
the quantum numbers S5 = C^^Qsi^ and = ip^C^/^Qsip* ■ Here Pa are traceless 
and Hermitian Nf x Nf matrices. Pa = p1 (a = 1, 2, • • • , Nj — 1), and Qs are symmetric 
Nf X Nf matrices, = Q5 (5 = 1, 2, • • • , A''/(A''/ + l)/2), in the flavor space. The chiral 
perturbation theory describing these NG modes for small ^ < mp/2 is exactly the same as 
that of SU(A''c) gauge theory with adjoint fermions at flnite fiB considered in because 
their symmetry breaking patterns are the same: the low-energy physics is dictated by the 
Dyson index /3 = 4. 

Let us consider the phase diagram of this theory. For m^/2 < /i < mp/2 at T = 0, it 
is energetically favorable for the U(l)s charged NG modes S5 with the excitation energy 
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— 2/i < to form the BEC, (S5) 7^ 0. On the other hand, at sufficiently large /i, 
the one-gluon exchange interaction in the ipip-channel is attractive in the color symmetric 
channel, leading to the condensation of the diquark pairing. Due to the Pauli principle, 
the BCS diquark pairing must be flavor symmetric, and takes the form {ip'^C^^Qs'tp) 7^ 0. 
This BCS pairing has the same quantum numbers and breaks the same symmetry as the 
BEC at small hb, and it is plausible that no phase transition occurs between the BEC and 
BCS regions. The phase diagram of this theory is similar to that of QCD at finite fii, as 
shown in Fig. |^. 

2.4 Sp(2A'^c) gauge theory at finite fiB = 1) 

We turn to Sp(2A'^c) gauge theory at finite ^b- From the property JcAf^Jc = (^4^^)* with 

Jc = -i(72(S)lN, (2.10) 



one has the relation: 



JcC^5'DCj5Jc = -V*, or [V,iJ,C-f5K] = 0. 



(2.11) 



The Dirac operator can be written in the form of real matrix [see ( [4. 4]) ], and the Dyson 
index is /3 = 1, the same universality class as two-color QCD Q ^] (see also ^^). One 
then finds that det[P(/i) + m] > for even Nf. Note here that even Nf is necessary for 
the positivity unlike S0(2A''c) gauge theory, because the quartet structure of eigenvalues 
(zA„, — iA„, zA* , — iA* ) reduces to one set of doublet {iXn, —iXn) when A^ is real or pure 
imaginary [^]. 

When m = hb = 0, because of the anti-unitary symmetry (2.11), chiral symmetry of 
the theory is enhanced to SU{2Nf) p8| , 49|. This can be seen by rewriting the fermionic 
part of the Lagrangian ( p.l| ) using the new variable ^ = {iPl^o'2T2'^*b)^ '■ 



(2.12) 
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7t 

{yj^C-f^JcQAij) = U(1)e restored 




which is manifestly invariant under S\]{2Nf). The chiral symmetry is spontaneously broken 
down to Sp(2A'^^) by the formation of the chiral condensate {"tpi^), giving rise to 2Nj — Nf — l 
NG bosons parametrized by the coset space S\J{2Nf)/Sp{2Nf). There are both \J{1)b 
neutral NG modes with the quantum numbers Ha = ipJ^Pa^P and V{1)b charged NG 
modes with the quantum numbers Y^a = ip^CjbJcQA'^ and = iP'^C'-^^JcQa^'* ■ Here Pa 
are traceless and Hermitian Nf x Nf matrices, 1,2,- •• ,A^/-1), andQ^ are 

antisymmetric Nf x Nf matrices, Q\ = —Qa = 2, • • • , Nf{Nf — 1) /2), in the flavor 
space. The chiral perturbation theory for small n < mp/2 is the same as that of two-color 
QCD at finite hb i, § because of the same symmetry breaking pattern dictated by the 
Dyson index /3 = 1. 

For m7r/2 < /i < mp/2 at T = 0, the U(1)b charged NG modes form the BEG, 
{Tia) 7^ 0. At sufficiently large fi, the one-gluon exchange interaction in the ^/''(/^-channel 
is attractive in the color antisymmetric channel and induces the condensation of the color 
and flavor antisymmetric BGS diquark pairing {ijj^ C^^JcQa^^) 7^ 0. Note that the diquark 
condensate has the different quantum numbers from the ones in S0(2A'^c) gauge theory. The 
BEC-BCS crossover of the diquark pairing is expected in this theory again, as depicted in 
Fig. I 

2.5 Brief summary 

Before proceeding, we summarize the results in this section. 

1. QCD at finite hb (or ///), S0(2A'^c) gauge theory at finite //b, and Sp(2A''c) gauge 
theory at finite ijlb belong to the different universality classes denoted by the Dyson 
indices /3 = 2, /3 = 4, and /3 = 1, respectively. 

2. QGD at finite pLj and S0(2A^c) and Sp(2A^c) gauge theories at finite ^b have no 
fermion sign problem and exhibit the BEG-BGS crossover phenomena in the phase 
diagrams. 
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Apparently, the phase diagrams of the three theories, QCD at finite fij, S0{2Nc) 
gauge theory at finite fiB, and Sp(2A''c) gauge theory at finite jiB, resemble each other 
qualitatively. In the next section, we will argue that these phase diagrams should be 
completely identical in the large- A'^c limit, including the chiral phase transition (not shown 
in figures). Also we will show that the phase diagram of QCD at finite /x^ should also be 
identical to them outside the BEC-BCS crossover regions. 

3. Orbifold equivalence in the large-Ai'c QCD and QCD-like theories 

In this section, we first briefly recapitulate the basic idea of the orbifold equivalence in 
large- Ac gauge theories. We then construct orbifold projection from S0(2Ac) or Sp(2Ac) 
gauge theory at finite /is to SU(Ac) QCD at finite fiB or fij. The relations between these 
theories via orbifold projections are summarized in Fig. 

3.1 Perturbative orbifold equivalence in pure gauge theories 

The main idea of the orbifold projection is as follows: we identify a discrete subgroup of 
the symmetry group of the "parent" theory, that is the S0(2Ac) or Sp(2Ac) gauge theory 
in our case. We then eliminate all of the degrees of freedom in the parent theory which are 
not invariant under the discrete symmetry. This gives the "daughter" theory, which will 
turn out to be the SU(Ac) gauge theory. We use a Z2 subgroup of the original S0(2Ac) 
or Sp(2Ac) gauge theory for the orbifold projection. In the large-Ac limit for fixed Nj 
(the 't Hooft limit), correlation functions of operators O^^^ in the parent theory invariant 
under the projection symmetry (which we call the "neutral operators"), and those of the 
operators O^^^ in the daughter theory made up of the projected fields, coincide to all orders 
in the perturbation theory |27| , 

{O^f^O^i^ ■■■), = {0?0i'^ ■■■),. (3.1) 

In this sense, the parent theory and the daughter theory are equivalent. Especially, mag- 
nitudes of neutral order parameters characterizing the phases should be exactly the same. 

Let us first consider the orbifold projection of the pure S0(2Ac) gauge theory. (For 
an earlier work of the orbifold projection from S0(2Ac) to SU(Ac), see |^].) For the 
"projection," we use Z4 subgroup of S0(2Ac) generated by Jc defined in ( p.lO| ). Jc satisfies 
the condition, 

tr( Jc") = 0, (3.2) 

when J" does not belong to the center of S0(2Ac), i.e., J" / ±l2Ar^. This condition is 
called the "regularity condition," which is utilized in the perturbative proof of the orbifold 

We can also construct the orbifold projections (not shown in Fig. |l|) from SO(2iVc) gauge theory with 
adjoint fermions at finite fiB to S\J{Nc) QCD with adjoint fermions at finite /is and to SU(A^c) QCD 
with antisymmetric fermions at finite fij, generalizing the argument of to the case with finite chemical 
potential. The resultant equivalence between SU(A'^c) QCD with adjoint fermions and SU(A'^c) QCD with 
antisymmetric fermions is a generalization of the onentifold equivalence |feo| to finite chemical potential. 
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equivalence. The transformation of induced by Jc is written as 

{Al%, ^ (Jc)n'(4°).','(J,)-], (3.3) 

which constitutes Z2 subgroup of the gauge group. We define the projection condition for 
the gauge field j4^^ to be invariant under the Z2 subgroup: 

{Al%j = {J,),,{Al%'r(Jc)j,]. (3.4) 

The gauge field jd^*-* satisfying this condition can be obtained by using the projector V 
defined as 

^(4°) = I E JcAfj- = \ {Af + JcAfj^') , (3.5) 

n=0 

where 4 is the number of elements of Z4. 

Remembering the property {A^'~')* = — the gauge field A^^ can be written as 

4S0 _ J A^ + - 
' ~ [C^ + D^^A^-B^)^ ^^-'^ 

where the fields A"^, B^, and (D^) are Nc x Nc anti-symmetric (symmetric) matrices. 
Under the Z2 symmetry, A^, are even while B^, are odd, so the orbifold projection 
sets B^ = = 0. Hence we have 

( zlA _T)S\ 

(3.7) 

If one performs a unitary transformation in the color space using the matrix 

Pc = ^( \ , (3.8) 

Afil^"^ can be expressed by the V{Nc) gauge field A^ = + ij4^ as 



(3.9) 



The top left component is the charge conjugation of the bottom right component, — {A^)'^ = 
(A^)'-^ . At large Nc, we can neglect the difference between U(A'^c) and SU(A''c) up to l/N^ 
corrections, A^ ~ A^^ Therefore, the action of the S0(2A'^c) gauge theory, 

Cso = -r^tr{F^?f, (3.10) 



^Notice the chemical potential is introduced as a boundary condition for the U(l) part of the gauge 
field at infinity. However it is not easy to impose an appropriate boundary condition at infinity in lattice 
simulations; although the integration region of the path integral must be limited so that it does not alter 
the boundary condition (in particular constant shift of the gauge field, — ^ + C, is forbidden), all 
field configurations, including constant shift, are summed over in actual lattice simulations. Because of this 
lattice artifact, /^-dependence disappears in \J{Nc) gauge theory |5^. In S\J{Nc) gauge theory, this problem 
does not arise. 
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is projected to that of the SU(A''c) gauge theory, 



^P-J = T^Mi^'-'')', (3.11) 

where F^^ is the field strength of the SU(A''c) gauge field ^4^^. In this way, SU(A''c) gauge 
theory is obtained from S0(2A'^c) gauge theory via the orhifold projection. 

We can similarly define the orbifold projection of the Sy>{'2Nc) gauge theory. The 
symplectic algebra Sp(2A'c) formed by 2Nc x 2Nc Hermitian matrices satisfies 



and can be written as 




^Sp ^ . (3.13) 



Here the fields {B^, C^, and D^) are Nc x Nc anti-symmetric (symmetric) matrices. 
If we choose the same projection condition for j4^p as ( |3.4D , 

Af = J,Afj-\ (3.14) 

one obtains B^ = = 0. This gives the SU(A'^c) gauge theory again. 

The mapping rule of the large- A'^c orbifold equivalence is as follows: if we equate the 
parent action with twice the daughter action for the Z2 orbifold projection |31, ^ ], 

'^so(Sp) 2£su, (3.15) 



these two theories are equivalent in the sense that ( p.lD holds within the neutral sectors 
i.e., the correlation functions of the neutral operators. The reason of this factor 2 will be 
clarified below [see ( [j.2lD ]. From the recipe ( [3.151) , we must take the action of the daughter 
SU(A^c) gauge theory as 



£su = j^tr(FSy)2, (3_ig) 



with the coupling constant gs\j satisfying 

530(3?) = 9S\J- (3.17) 

Let us give a schematic of the proof of the orbifold equivalence to all orders in the 
perturbation theory following Bershadsky and Johansen . This can be easily generalized 



to the cases with the fundamental fermions [35|, as we will see in the next subsection. We 



assume that the gauge fixing condition is consistent with the projection - the ghosts are 
related by the projection, and all propagators in two theory takes the same form up to the 
color factors. 

As a pedagogical demonstration, consider a vacuum planar diagram of the S0(2A'^c) or 
Sp(2A^c) gauge theory in Fig. ^. In order to obtain the SU(A'^r;) diagram from the S0(2A^c) 
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or Sp(2A^c) diagram, we insert the projector V to each propagator. Because the 't Hooft 
coupUngs are taken to be the same and the propagators are the same up to color factors, 
the only difference, if exists, comes from the contractions of color indices. Remembering 
( |3.5| ), the additional kinematic factor multiplied by the SU(A'^r;) diagram is 

/I \ A'p 

^ - • tr(J-"i J"5) • tr(J-'^V-"V6) • tr(J-"V-"V-"«) • tr( J"i J"V^»), 

n.=0,l 

(3.18) 

where J = —ia2 is a 2 x 2 matrix and the factor (1/2)^^ comes from the projectors with 
Np = 6 being the number of propagators. Because of the regularity condition ( |3.2D , it is 
nonvanishing only when 

(3.19) 

The last condition follows from the others, and hence, there are A'^^^ — 1 = 3 independent 
constraints, where A'^^^ = 4 is the number of color index loops. In ( |3.18| ) the combinations 
of nj = 0, 1 for all Np = 6 propagators under the Ni — 1 = 3 constraints give a factor of 
2^~^, traces over color indices give 2^, and thus, the total factor is 

2-6 . 2^-3 . 2*^ = 2. (3.20) 

Generally, for any given planar vacuum diagram with Np propagators and A'^^ loops, 
the projectors give a factor of (1/2)'^^, the combinations of nj = 0, 1 (i = 1, 2, • • • ,Np) 
under the N^ — 1 constraints give 2^^^ and the traces over color indices give 

Therefore, the total factor is always the same: 

2-7Vp . 27Vp-(7Vi-i) . ^Nl = 2. (3.21) 

This factor 2 reflects the fact that the number of degrees of freedom in the parent theory 
is twice larger than that in the daughter theory. Hence the vacuum energy per degree of 
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freedom is identical between these theories. This is why the recipe ( 3.15| ) is necessary to 



match the degrees of freedom between the parent and daughter theories. However, this 
argument does not hold for nonplanar diagrams. Actually, one can check that the number 
of independent constraints is no longer — 1, and the kinematic factor is different from 
2 counted in ( 3.21| ) p7| . This is why we need to take the large- A'^c limit to suppress the 



nonplanar diagrams. We can repeat the same argument to any correlation functions of 
neutral operators (gauge-invariant operators in this case) in any dimension and any gauge 
group. 

Note that this equivalence holds as long as the projection symmetry, Z2 subgroup of 



S0(2A^c) X U(1)b or Sp(2A^c) x U(1)_b, is not broken spontaneously [Q. We will come to 
this issue in more detail in Sec. 13. 



3.2 Perturbative orbifold equivalence with fermions 

In this subsection, we further include the fermions and define the orbifold projections from 
S0(2A'c) and Sp(2A''c) gauge theories to SU(iVc) gauge theory at finite or ^j}^ 

In order to obtain fermions at finite /is, we use Z4 subgroup of S0(2A'c) or Sp(2A'c) 
gauge group generated by Jc and Z4 subgroup of U(1)_b generated by w = e*'^/^. We choose 
the projection condition as 

V.S0(SP)=^(JJ,,,V.'P^SP), (3.22) 

which generates Z2 subgroup of S0(2iVc) x U(1)b or Sp(2iVc) x U(l)ij. The color 2Nc- 
component fundamental fermion is decomposed into two A^c-component fields, 

with a and b being the color indices. Performing the unitary transformation Pc defined in 
( |3.8| ), we have 

p^^so(Sp) = ( ^+ ) , (3.24) 




where il}± = {ipa i iipi,)/\/2. The fermions ip+ aiid which couple to (AJ^^)^ and 
from (|3.9| ), transform as antifundamental and fundamental representations under SU(A'^c); 



respectively. After the projection ( 3.22 ), only the fermion ip-. survives. Taking into account 



the relation ( 3.15 ), the action of the daughter theory is given by 



^su = j\- tr(FSU)2 + J2 V^f {I'-D, + m + /.7^) V^f , (3.25) 



^''The orbifold projection can be generalized to the case with both finite jj.B and finite fij. More generally, 
one can construct the orbifold projection in the system where each flavor has different chemical potential 
M/ if = 1,2,- ■■ ,Nf). 



-15- 



where -00^ = V^V*-) = + iA^ . This theory is QCD at finite baryon chemical 
potential //_b = A^cA*- 

On the other hand, in order to obtain fermions at finite fij for even Nj, we use 
Jc G S0{2Nc) [or Jc G Sp(2A^c)] and Ji £ SU(2)isospm defined by 

Ji = -ia2 (8) lNf/2- (3.26) 
We choose the projection condition to be 

The flavor A^j-component fundamental fermion is decomposed into two (A^j/2)-component 
fields, 

V.SO(Sp) = (^^ ^^)^ (3.28) 
with / and g being the flavor indices. If we define (f± = =F ^V'±)/V2 and .^-i- = 



('i/'i ± i'0±)/V2, the fermions (/3± survive but ^± disappear after the projection ( 3.27] ). 
Because ip± couple to (A^)'-' and A^ respectively, the action of the daughter theory is 
expressed as 

Csv = -\- tr(FSU)2 + ^ ^su + ^ i ^su^ (3 29) 

where = V^f- and -0-^ = V^^+- This theory, in which "zb" can be regarded as the 
isospin indices, has the isospin chemical potential fij = 2/u. 

Given the orbifold projections, we now prove the orbifold equivalence of the gauge 
theories with fermions at finite chemical potential. Let us consider a diagram (a correlation 
function or an observable) in S0(2A'^c) or Sp(2A'c) gauge theory which have quark loop(s). 
Here we write both color and flavor index lines for each quark line since the contractions 
of flavor index loops also give kinematical factors. 

We first look into the case with finite /i/. Because the matrix Jc multiplied from the 
right in ( \3.4i ) is just replaced by Ji in ( 3.27| ), we can repeat the proof for the pure gauge 



theory straightforwardly. For a diagram with Np propagators and A'^'^^ color (A'^^'^^ flavor) 
index loops, we count additional kinematical factors originating from contractions of both 
colors and flavors that are multiplied by the parent SU(A'c) diagram. The projectors for 
the propagators give a factor of (1/2)^^, the combinations of rij = 0, 1 (i = 1, 2, • • • , Np) 
under the A'^^'^^ + A'^^'^'' — 1 constraints give 2^^~(^i ^+N^^^-i) ^ ^^^^ ^j^g traces over color and 
flavor indices give ' and 2^l^ respectively. Therefore, the total factor is 

. . 2< . 2^i'' = 2, (3.30) 

( f) 

for any number of quark loops A'^^ . 

On the other hand, this argument is not applicable to the case with finite ^p since 
the structure of the projection condition ( |3.22| ) is different from (|3.4| ) and ( 3.27] ). However 



we can still justify the orbifold equivalence as long as the number of quark loop is one. 
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(c) 

Consider a diagram with Np propagators and A'^^ color index loop and one fermion index 
loop. The fermion index loop gives no condition for J"'' in Eq. (3.19), and the number 

(c) 

of the constraints for Tij = 0, 1 (i = 1, 2, • • • , Np) is Nj^ . Then one again concludes that 
the total kinematical factor is 2. Note that this does not hold any more if the number of 
fermion index loop is more than one. 

Therefore, remembering that the fermion loops are suppressed by a factor of Nj/Nc, 
the orbifold equivalence at finite fi b holds only in the 't Hooft limit (large- A^c limit for fixed 
Nf) while that at finite /// holds both in the 't Hooft and the Veneziano limits (large- A^c 
limit for fixed Nf/Nc). 

One can also argue to what order of 1/Nc the orbifold equivalence is valid. The leading 
corrections to the 't Hooft limit come from the one-fermion-loop planar diagrams, which, 
as we have shown above, do not distinguish between jjLp and /i/. Therefore, the difference 
of the expectation values of gluonic operators, such as the Polyakov loop, between QCD 
at finite ^p and /i/ is at most ~ (Nf/Nc)'^ due to a two-fermion-loop planar diagram. In 
particular, the difference of the critical temperature of the deconfinement is ~ {Nf/Nc)'^, 
as previously discussed in [52] without using the orbifold equivalence. 



3.3 Condition for the orbifold equivalence and the BEC-BCS crossover region 

From the orbifold projections we have constructed, the expectation values of the neutral 
operators for given T and fi should be identical between the parent and daughter theories. 
Note here that, in order for the orbifold equivalence to QCD at finite fip to hold, the 
projection symmetry, Z2 subgroup of the S0(2A^c) x U(1)b for the S0(2A^c) gauge theory, 
or Z2 subgroup of the Sp(2A'c) x \J{1)b for the Sp(2iVc) gauge theory, must not be broken 
spontaneously ||3^. However, the U(1)b symmetry, whose Z4 subgroup is used for the 
projection to QCD at finite jip, is broken down to Z2 inside the BEC-BCS crossover 
regions in figures |3| and ^ Therefore, the orbifold equivalence to QCD at finite ^p is valid 
only outside the BEC-BCS crossover regions. Note that, even at large //, above the critical 
temperature of the superfluid, T > Tc, the condensates melt away and the equivalence 
recovers. 

On the other hand, the equivalence to QCD at finite m is valid as long as the projection 
symmetry, Z2 subgroup of S0(2A^c) x SU(2)i;' or Z2 subgroup of Sp(27Vc) x SU(2)ir, is not 
broken. This condition is satisfied for arbitrary // and T in all the three theories with flavor 
symmetry, QCD at flnite fij and S0(2A''c) and Sp(2A''c) gauge theories at finite fip. As a 
result, the chiral condensate {^pijj), the superfiuid gap A, and their critical temperatures are 
identical; the phase diagrams characterized by these order parameters completely coincide 
in the large- A^c limit. 

In particular, QCD at finite fip is equivalent to QCD at finite fj-i outside the BEC- 
BCS crossover region. Remembering that dropping the complex phase of the fermion 
determinant of QCD at finite fj,p reduces to QCD at finite /U/ for Nf = 2, this means that 
the phase- quenched approximation in this region is exact at large-Nc for neutral operators, 
e.g., the chiral condensate. This explains why the phase-quenched approximation works in 
QCD at high temperature, as already observed in the lattice QCD simulations for small 
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At 



and in model calculations such as the chiral random matrix model |20|, Nambu- 



Jona-Lasinio model [21, and hadron resonance gas model |22]. As we will show in 



sections 4.2 and 4.S. the results of the chiral random matrix model can also be understood 



as a consequence of the orbifold equivalence in this model. 

3.4 Nonperturbative orbifold equivalence at high density 



In Sec. 3.1, we discussed the orbifold equivalence at the perturbative level in the large 
Nc limit. One might ask if the equivalence does hold nonperturbatively and to what 
extent the 1/Nc corrections are important. For SU(A''c) QCD at large m and S0(2A'c) and 
Sp(2A^c) gauge theories at large fiB, '^^ can answer both questions for several quantities 
explicitly. This is because the coupling constants ^sui 5sOi and i^sp are small due to 
the asymptotic freedom and the calculations are under theoretical control. Note that the 
weak-coupling does not necessarily mean that calculations are perturbative: as shown in 
( 3.35 ), ( p. 36 ), ( p. 37 ), and ( 3.3^ ) below, the coupling constant dependences of the BCS 
gaps are indeed nonperturbative. Although our calculations do not constitute the proof 
of nonperturbative orbifold equivalence for all the neutral observables, they provide a 
nontrivial piece of evidence for the equivalence. 

Let us consider Nj = 2. First of all, the chiral condensate {i/^ip) is vanishing in all the 
theories at asymptotically large /i, and the equivalence is trivially satisfied. The nontrivial 
quantity we can compare is the superfluid gap A. Its equivalence is required, e.g., from 
the equivalence of the fermion occupation number ('^^^(t, x)'i/;bg(t, y)). The gap A can 
be computed using the technique in [p^. The main modification in the gap equations 



compared with [53| is the group theoretical factor involving Nc given by 

iVc + 1 



/rpSV\T (rpSlJ\ |'T^SU\ 



2Nr. 



■SU^ 



\2Nc 
13-y y-'-a )lS 



SO^ 



2A^c 

2Nr 



-{tT. 



aS 1 



1 



2Nc 
5 



4 

2Nc + l 



(3.31) 
(3.32) 
(3.33) 
(3.34) 



respectively. Solving the gap equations, we obtain the following BCS gap (up to prefactor) 
for each theory: 



~ /"exp ( --— 



A^u ^ 

Mi 



//exp 



~ /"exp 
A^P ~ ^exp 




(3.35) 

(3.36) 

(3.37) 
(3.38) 
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Equation ( 3.35 ) is the result obtained in [^4|. For Nc = 3, and reduce to the 



results obtained in |Q and respectively. Note that in the 't Hooft limit A^y, A^^, and 
A^^ remain finite (and the BEC-BCS crossover regions of these theories do not disappear) 
while A^y is vanishing. This originates from the fact that the diagrams of the one-gluon 
exchange responsible for A^y, A^^, and A^^ are planar, whereas it is nonplanar for A^^ 
| [55| , p6| . These consequences are consistent with our claim in Sec. that the orbifold 
equivalence holds between QCD at large fij and S0(2A''c) and Sp(2A^c) gauge theories at 
large fiB while the equivalence with QCD at large fiB does not hold inside the BEC-BCS 
crossover region. This gives a simple example in QCD that, if the projection symmetry is 
broken, the orbifold equivalence is not valid. 

Now let us compare A^y, A^^, and A^^. For the comparisons to the leading order, 
we take the ratios of the factors in the exponential between ( 3.36| ), (3.37), and ( 3.3^ ). 
Remembering gs\j = gso = dSp^ the ratios read 



( N ^ 2(iV2 - 1) I 0.873 (iV, = 3) 

asp/su(A^c) = ^^^(2Ar, + l) = \ 1 (iV, = oo) " ^'-'^^ 

Clearly, the equivalence holds in the large- A^c limit between QCD at large fij and S0(2iVc) 
and Sp(2A^c) gauge theories at large fiB- It also turns out that the orbifold equivalence 
is rather well satisfied even in real QCD, Nc = 3. For the complete equivalence of the 
BCS gap, one has to check if the prefactors in (|3.36|) , (|3.37D , and (|3.38| ) are also identical, 
for which subleading effects are important. Here we simply assume the equivalence of the 
prefactor in the large- A'^c limit, and see whether the equivalence of other nonperturbative 
quantities follow or not. 

Provided A^^/A^y — )- 1 and A^^/A^y — )- 1 in the large- iVc limit, the critical temper- 
atures Tc of the superfluid phases also coincide with each other, since Tc is proportional to 
the BCS gap with the universal proportionality factor: 

Tc = -A, (3.41) 
vr 

where 7 ~ 0.577 is the Euler-Mascheroni constant. 

Moreover we compute the magnitudes of the diquark and pion condensates following 



|5j] as 



/J^, ,,\SU jSU ijSUi o / 6A^c M^A^y /lo^ 

(d75n),, = d,, , Id,, I = 2^ -^2^^^, (3.42) 

ii'lC.^i.lfO = s^,8fHl% \dl-\ = 2/^^^, (3.43) 

{¥aCl,rA = {Jc)aMy^d%, \d%\ = 2^^^^^, (3.44) 
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where Jc and Jj are defined in ( |2.10D and ( 3.26 ). The ratios read 



l^syi - «so/su(iVc)^, (3.45) 

I jSp I A Sp 

1^ = «Sp/su(iVc)^, (3.46) 



where aso/su(-^c) and a^^i^ulN,.) are the quantities defined in ( p. 39 ) and ( |3.40| ). Both 



indeed approach unity in the large- A'c hmit when A^^/A^y — )• 1 and A^^/A^y — t- 1. 

Although we did not attempt in this paper, one should also be able to check the 
equivalence for other quantities at large ;U, such as the four-quark condensate {{il)tp){'4)tp)) , 
the pion decay constant /tt, and so on. 

3.5 Brief summary 

We summarize our results in this section. 



1. The whole phase diagrams described by the chiral condensate and the superfluid gap 
should be universal in the large- A'^c limit between QCD at finite /i/ and SO (2 Ac) and 
Sp(2Ac) gauge theories at finite ^b- 

2. The phase diagram of QCD at finite should also be identical to those of other 
theories above outside the BEC-BCS crossover regions. In particular, the phase- 
quenched approximation for the chiral condensate is exact in this region. 

3. At asymptotically large chemical potentials, the equivalence is rather well satisfied 
even for Ac = 3. From this fact, the phase-quenched approximation for the chiral 
condensate is expected to work well in real QCD. 



4. Orbifold equivalence in the chiral random matrix theories 

If the orbifold equivalence holds between the original gauge theories, it is natural to expect 
that the equivalence should hold at the level of the corresponding low-energy effective 
theories. In this section, we show that the orbifold equivalence holds perturbatively as 
well as nonperturbatively in the chiral random matrix theory (RMT), a solvable effective 
theory of QCD (or QCD-like theory) first introduced in [^7|. Actually the RMT has the 
size of the matrix A, which is taken to infinity (thermodynamic limit) in the end. In 
this sense, the RMT is a "large-A" matrix model, and hence, the perturbative proof of 



the orbifold equivalence given in Sec. 3.1 is applicable. Note that the size of the random 



matrix is not related to the number of color Ac. The orbifold equivalence will be verified 



nonperturbatively in Sec ^ by solving the RMT following |20|, |58|. 
4.1 Chiral random matrix theories 

In this subsection, we briefly review the basic aspects of the RMT. For reviews of the RMT 
in more detail, see Isl, |39| . 
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4.1.1 Chiral random matrix theories at small chemical potential 

The partition function of the RMT is given by an integral over a Gaussian random matrix 
ensemble, instead of the average over the gauge field of the original Yang-Mills action: 

Z = I dcDjJdetPe-^^'*^-*'*, (4.1) 

i=i 

where $ is an x [N + u) random matrix element, is the size of the system, and v is 
the topological charge. We also introduced a suitable normalization with the parameter 
G in the Gaussian. Note that there is no spacetime coordinate in the theory; the size 
of the matrix N corresponds to the spacetime volume. It is taken to infinity in the end, 
corresponding to the thermodynamic limit. 

The matrix structure of the Dirac operator T) is chosen such that it has the correct 
anti-unitary symmetries and it reproduces the correct global symmetry breaking pattern 
of the system. In particular, anti-Hermiticity and chiral symmetry of the Dirac operator 
at fi = m = require 

pt = _p, {V,j^} = 0. (4.2) 



The quark mass m, quark chemical potential /U |5^, and temperature T can be incor- 
porated into D. While T does not destroy any relation in ( |4.2| ), /x and m break the former 
and the latter relations in (|4.2| ), respectively. 

Depending on the anti- unitary symmetries of the Dirac operator, the ensemble is dis- 
tinguished with the real, complex, or quaternion real [see ( [4.14 ) for the definition] matrices 



denoted by the Dyson index /3 = 1, /3 = 2, and /3 = 4, respectively. The value of /3 cor- 
responds to the degrees of freedom per matrix element. QCD or QCD-like theory in each 
universality class and the corresponding RMT is listed as follows pT| : 

• The Dirac operator of SU(A'^c > 3) QCD has no anti-unitary symmetry and the 
corresponding RMT belongs to /3 = 2. The Dirac operator is taken as 

V=( 7' ^ + ^'V (4.3) 
\^-«>t + ^l m/1 J' ^ ^ 

where ^> is an x (N + u) complex matrix and mj (/ = 1, 2, • • • , Nj) are the quark 
masses. 

• The Dirac operators of SU(A''c = 2) QCD and Sp(2A''c) gauge theory have the anti- 
unitary symmetries, PT2^75C = ^^CV* and ViJc^^C = iJc^i^CV* , respectively. 
Here C is the charge conjugation matrix, ^ is the antisymmetric generator of SU(2), 
and Jc is defined in ( 2.10| ). The corresponding RMT belongs to the universality class 



/3 = 1 and the Dirac operator is taken as 



m/1 <1> + /xl 
where <I> is an x (A^ -|- v) real matrix. 



V= 7 , (4.4) 

l-$^ + /il m/1 j ^ ' 
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• The Dirac operators of SU(iVc) QCD with adjoint fermions and SO(A'^r;) gauge theory 
have the anti- unitary symmetry, P75C = 75 CP*. The corresponding RMT belongs 
to the universahty class /3 = 4 and the Dirac operator is taken as 

V=( 7' ^ + ^'V (4.5) 

where $ is a 2N x 2(A^ + u) quaternion real matrix [see ( [4.14] ) for the definition]. 
These RMTs at finite fi can alternatively be formulated by the two-matrix representation 



62 , 63 1 where the identity matrix multiplied by // is replaced by the random matrix element. 
The effect of temperature T can be incorporated as the (first) Matsubara frequencies by 
changing ^ ^ 11 + iT for one half of the determinant and ^ ^ ^ — iT for the other half of 
the determinant in the simplest model |60||.^'^ 

There is a regime (called the e-regime) where the RMT is exactly equivalent to QCD: 
when the typical scale of the system L is much smaller than the pion Compton wavelength 
and is much larger than the inverse of the scale of chiral symmetry breaking pi 



< L < — , nL < 1, (4.6) 

QCD reduces to a theory of zero momentum modes of pions. In this regime, the system 
has a universality, i.e., the dynamics depends only on the symmetry breaking pattern and 
is independent of the microscopic details; QCD can be replaced by the RMT with the same 
global symmetry breaking pattern. Outside the e-regime, the universality is lost. However, 
the RMT is still useful as a schematic model to study the qualitative properties of QCD 
such as the phase structure at finite T and /i |^^. The advantage of the RMT is that it 
can be solved analytically although QCD itself cannot be. 

4.1.2 Chiral random matrix theories at large chemical potential 

Recently, a new class of RMTs which describe the superfluid phase at asymptotically large 



// have been identified in [36|. The partition function is given by 

Nf 

Z = I dAdBl[detV ^-^GHHA^A+B^B)^ (47) 

i=l 

where A and B are N x N spacetime independent random matrix elements. Here only 
the topological sector 1/ = is considered because the topological susceptibility is strongly 



suppressed at large fi [69, 70 [. Chiral symmetry 

m75} = 0, (4.8) 



^^Introducing T in this way may break fiavor symmetry of the system in some cases. However, the final 
result is shown to be equivalent to that of the correct prescription preserving flavor symmetry j6^] . 
^^The mathematical aspects of the same two-matrix model for /? = 1 were previously studied in Js?! , 



- 22 - 



is preserved at finite // in the chiral limit m = 0, but anti-Hermiticity of tlie Dirac operator 
= -V is lost. The non-Hermitian Dirac operator is taken as 

V^H'" ), (4.9) 

where both A and B are the real, complex, or quaternion real matrices, denoted by the 
Dyson index (3 = 1, (3 = 2, and /? = 4, respectively. The /? = 1 RMT corresponds to SU(2) 
QCD H] and Sp(iVc) gauge theory at large hb, (3 = 2 RMT to SU(A^c > 3) QCD at large 
m, and (3 = 4: RMT to adjoint QCD and S0(2A''c) gauge theory at large ^b- 

We can define the e-regime at large fi where the system has the universality and QCD 



is equivalent to the RMT |7|]: 



^«L«— , (4.10) 

where A is the BCS gap and ttZtt is the pion mass associated with chiral symmetry breaking 
by the diquark condensate (not by the usual chiral condensate) . It is shown in [^6| that the 
partition function of the RMT actually coincides with that of the finite-volume effective 
theory of QCD at large fi. 

4.2 Orbifold projections in the chiral random matrix theories 

In this section, we construct the orbifold projections in the chiral random matrix theories 
(RMTs) between /3 = 4, (3 = 2, and (3 = 1. Thereby a class of observables in the RMTs 
with the different Dyson indices are found to be identical to each other. In the following, 
we will concentrate on the RMT at finite // and T = introduced in Sec. |4.1.1| , which can 
be easily generalized to nonzero T. For simplicity, we set u = and consider degenerate 
quark masses mj = m. The generalizations to the high-density RMTs in Sec. 4.1.2| and to 



the RMTs in the two-matrix representation [S2, 63] are straightforward. 



The construction of the orbifold projections is as follows: we start with the /3 = 4 or 
(3 = 1 RMT at finite hb with the size of the matrix $ being 2A^, and define the orbifold 
projection to the (3 = 2 RMT at finite hb or /i/ with the size of the matrix A^. The relations 
between these RMTs via orbifold projections are summarized in Fig. ^.^^ Each orbifold 
projection reduces the independent degrees of freedom of the theory to half. 

4.2.1 Orbifold projection from (3 = A to (3 = 2 

The partition function of the (3 = A RMT is given by 

Z = I dM^ e"^, S = Sb + Sf, (4.11) 



where 



Nf 



Although not shown explicitly in this paper, one can also construct the orbifold projection from the 
13 ^2 RMT at finite fj.i with the size iV to the /3 = 1 RMT at finite fis with the size iV. 
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j8=4 or j8=1 RMT with jJB 



Equivalent outside 
the BEC-BCS 
crossover region 



Equivalent 
everywhere 



/3=2 RMT with ijb 



► 



13=2 RMT with jji 



Equivalent outside the 
BEC-BCS crossover region 

Figure 6: Relations between /3 — 2 RMT at finite or /i/ and /3 = 4 and (3 = 1 RMTs at finite 
/is. (3 = 2 RMT at small and large /i/ can be obtained from (3 = A and (3 = 1 RMTs at small and 
large through orbifold projections. (3 = 2 RMT at small ^b can also be obtained from (3 = A 
and (3 = 1 RMTs at small ^b while (3 = 2 RMT at large /is (inside the BEC-BCS crossover region 
of the parent RMTs) cannot be. 



and 

V=( "^'^^ '^ + f^^2N\ (4^3) 
y-^^ + fll2N ml2N J 

Here are complex Grassmann 4A^-component vectors and <I> is a 2N x 2N quaternion 
real matrix of the form: 

where are N x N real matrices and = {—i,ak) with Pauli matrices a^- 
For the bosonic matrix <I>, we impose the projection condition 

Then we obtain 

(4.16, 

which is equivalent to two copies of a x complex matrix after a unitary transformation 

where cj) = a^+ia^- The bosonic part of the action is mapped to the one for the (3 = 2 RMT. 
Note that the factor 2 in the recipe ( 3.15| ) is reflected in the difference of normalization 
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between /3 = 4 and /3 = 2 in ( [4.12| ) if the trace for /3 = 4 is understood as the so-cahed 
"QTr" which is one-half the usual trace. 

In order to introduce a projection for the fermions, we write ^ by using two 2N- 
component fermions and V'L as 



Here V'ij and ipL are further decomposed into two A^-component fermions 



(4.18) 



(4.19) 



where the flavor index is suppressed for simplicity. The projection to the (3 = 2 RMT at 
finite //B is given by 

lpR = UjJ'lpR, tpL=i^JipL-, (4.20) 
where oj = e™/"^ as defined in Sec. 3.2. Performing the unitary transformation, we have 



'R 



^R 







(4.21) 



where = {ip}^ it iipj^) / y/2 and ■i/'^ = (V'i =t i'4)\)/V2. The fermionic part of the action 
reads 



R) 













mlN 





















\ / \ 


mlN / \i^L J 



where 



^(/^)/3= 



1 



13=2 



R 



\/2 I V'i 



2^'/3=2^?(A*)/3=2^'/3=2, 

(4.22) 

(4.23) 



This is the /3 = 2 RMT at finite quark chemical potential [i. 

In order to obtain the (3 = 2 RMT at finite /i/, we impose the projection condition: 

J^rJ;^ = ^R, JMr^ = ^L. (4.24) 

If we define (p^ = (V'^j =F ^V'^g)/v^ and = (iI'rj ± iil)^g)/^/2 (and similarly for (^9^ and 
^f^) with / and g being the flavor indices, the fermions survive but S^^^^ disappear 
after the projection. The fermionic part of the action reads 





6^ + fllN 





( 


i>* + /"liv 





\ 








ct 


) + /xlAr 




fR 














(j)^ + hIn 







ml]\f 


J 




^?(/«)/3=2^^ 


J=2 











(4.25) 
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where 

(4.26) 





Because and have the quark chemical potential +// and respectively, it is the 
/3 = 2 RMT at finite isospin chemical potential /U/ = 2fi. 

4.2.2 Orbifold projection from /3 = 1 to /3 = 2 

The /3 = 2 RMT can also be obtained from the /3 = 1 RMT. We start with the action of 
the /3 = 1 RMT given by ( 4.12| ), but the Dirac operator is now 



y-^^ + Hl2N 'rnl2N J 
where $ is a 2N x 2N real matrix. <I> can be parametrized as 

\ —a + a a — a I 

where are N x N real matrices. Note that the only change in this expression compared 
with ( 4.14 ) is that the factors i in front of and are absent. Then one can easily find 



that the same projection conditions for <I> and as Sec. [4.2. 1| give the (3 = 2 RMT at finite 
fi or finite /i/. 

4.3 Solving the chiral random matrix theories 

The orbifold equivalence in the RMT predicts that the /3 = 4 and /3 = 1 RMTs at finite fiB 
and /3 = 2 RMT at finite fij are equivalent to each other in the neutral sector. Moreover, 
outside the superfluid phase, the above three theories must also be equivalent to the /3 = 2 
RMT at finite ^b- In this section, we will verify these statements at the nonperturbative 
level by computing the effective potential of each RMT. 



For the j3 = 2 and /3 = 1 RMTs, the effective potentials are computed in pO| , 58]. The 
equivalence of the effective potential of the j3 = 2 RMT at finite ^ b (at = 0) and that of 
the j3 = 2 RMT at finite (at = 0) is pointed out outside the pion condensed phase. 
Here we show that the equivalence holds between a larger class of RMTs as a natural 
consequence of the orbifold projections constructed in the previous subsection. 

Let us sketch the derivation of the effective potential of the RMT. First, introduce 
the Grassmann vectors (fermions) tp to write the determinant into an exponential form. 
Second, perform the Gaussian integration over the matrix element which leads to the 
four-fermion term in the exponent. Third, introduce the bosonic auxiliary field A to make 
it the fermion bilinears (this procedure is called the Hubbard-Stratonovich transformation) . 
Fourth, perform the Gaussian integration over ^ . Finally, the effective potential is given 
as the saddle point of the integrand in the — )• oo limit (the thermodynamic limit). 



^^Our arguments depend on the ansatz of the saddle point of A (defined below) at finite T and /i similarly 
to @||. 
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In the following, we consider the Nf = 2 RMT in ( |4.l| ) with the quark mass m f and the 
chemical potential fif for each flavor, / = 1,2. The baryon and isospin chemical potentials 
are defined as 

= ^ = ^(/"i+M2), (4.29) 

= Y = ^(m-/^2). (4.30) 

We denote the chiral condensate as crj, pion condensate as p, diquark condensate as A, 
and their sources as mj, A, and J, respectively. 

4.3.1 Effective potential of /3 = 4 

We first consider the /3 = 4 RMT with degenerate quark mass mf = m at finite baryon 
chemical potential pf = fls (and hence o"/ = a). We will focus on T = and introduce T 
later changing as p-B — ?• fis + iT and pB fiB — iT for each half of the fermion determinant 
following 1 64]. 

Remembering the definition of the quaternion real matrix ( [f.l4 ), the partition function 
can be rewritten as 

Z = j da^'d^*d'^ exp [-2iVG2 tr(a^a^) + i^^RMij^a^ + flBhj)i^ij (4-31) 

where ipR^i are the Grassmann 2A^-component vectors and f,g are flavor indices and i,j 
run over 1, 2, • • • 2N. Integrating out and using the Fierz identity, 

(^m)"'«)"' = ^^acSdb, (4.32) 

the partition function reduces to 



/ 



dip* dip exp 



(4.33) 



Performing the Hubbard-Stratonovich transformation by introducing the auxiliary real and 
symmetric Nf x Nf matrices Kfg and Lfg, one has 



dKdLd^P*d^ exp [-8NG^ tv{K^ + L^) + fiBii^^^i^li + 

+'^^m^miK + ^L)fg + 2i:%i^{^{K - iL)fg 
+Ml^P{^A + Mfg4^{*rLj\ . (4.34) 
Integrating over the fermionic variables tp and ip* leads to the expression: 

Z = J dAexp[-Nnp=4{A,A^)], (4.35) 



-27- 



where ri/3=4 is an effective potential given by 

rj/3=4 = SG^ ti{AA^) - 2 log det Q. 
Here A = K — iL \s the complex and symmetric Nf x Nf matrix and 



Q 



flsSfg 2A + M 



(4.36) 



(4.37) 



We set the source term M and make the ansatz for A as follows, 

M 



1 m 




A = 










r 




\ iJ 


m 1 




[ iA 





(4.38) 



Shifting a and A such that m and J dependences are absorbed into the quadratic term 
and adding the T-dependence, one finally arrives at the effective potential of /3 = 4 RMT: 



/3=4 



a--] + A-- 



2^1n[4a2 + 4A2 - (/XB ±iT)2]. (4.39) 



The chiral condensate and the diquark condensate are expressed using a and A as 



{uu)p=4^ = — 9mlnZ^=4 



m=0 



J=0 



-4GV;3=4, 



-4G2a^=4. 



(4.40) 
(4.41) 



4.3.2 EfTective potential of /3 = 2 

Similarly to the case with /3 = 4, one can obtain the effective potential of /? = 2 RMT. 
This was previously computed in |2C] and the result reads 

n^=2 = G\ai - mif + {02 - m2f + 2{p - Xf] 

- i ^ ln[(ai + ^1 ± iT) ((72 -fi2T iT) + p^] [{ai -piT iT) {a2 + ^2 ± iT) + p^] . 

(4.42) 



The chiral condensate and pion condensate are calculated as 

1 



{uu)fi=2 = J^^nii lnZ;3=2 



mi=0 



{d'j^u)i3=2 = ■^dxlnZi3=2 



(4.43) 
(4.44) 



A=0 



Note that, as long as p = (i.e., outside the pion condensed phase), the potential (|4.42 ) is 
a function of p\ = {pB + pi)"^ and p\ = [ps — pi)"^ ■ This property leads to the relation: 



^^/3=2(/^b)Uj=0 = f^/3=2(/^7-)Ufl=0 for p = 0. 



(4.45) 
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Here p = is the condition that the projection symmetry, which is used for the orbifold- 
ing in Sec. 4.2, is not broken spontaneously. From ( [4.45 ), the magnitude of the chiral 
condensate a and the critical temperature of chiral phase transition T°' in each theory 
coincide, 



(7/3=2 (/"B ) U/ =0 = 0-/3=2(/^/) Us =0 for p = 0, 
Tp=2ifiB)\fli=0 = T^=2ifil)\flB=0 foi" P = 0, 



(4.46) 
(4.47) 



as a consequence of the orbifold equivalence. Especially, this shows that the phase-quenched 
approximation for a 13=2(^3) and T^^2iP'B) works outside the pion condensed phase, as 
mentioned in |20]. 

It should be remarked that, even though the effective potentials are identical in ( 4.45| ) 
for p = 0, the partition functions themselves are not generally the same. This is because 
the preexponential factor also contributes to the partition function, which is not taken into 
account in computing the effective potential. Therefore, the sign problem measured as 
the phase of the partition function can be severe inside as well as outside the pion condensed 
phase |^3|. The result here shows that the phase-quenched approximation is exact for the 
observables above independently of the severity of the sign problem, as long as p = 0. 

4.3.3 Effective potential of /3 = 1 

The effective potential of /3 = 1 RMT is computed in Kq] as 



/3=1 



G^[{ai - mif + (as - ms)^ + 2{p - A)^ + 2(A - jf 
1 



ln{[{ai + pi ± iT){a2 - fi2T iT) + + A^] 

X [(cTi - /ii ± iT){a2 + /i2 =F iT) + p^ + A^] + 4AVi/i2} 

x{[(ai - pi ^ iT){(j2 + P2 ± iT) + p2 + A^] 

X (0-1 + pi =F i'T){a2 - P2 ± iT) + p^ + A^] + 4A2pi^2}. 

The chiral condensate, pion condensate, and diquark condensate read 

1 



{uu)p=i 



2N 

1 
4iV 

1 

m 



dmi lnZ^=i 



:1) 



rrai=0 



A=0 



-G^A 



/3=1- 



J=0 



The potential ( 4.48| ) has the symmetry 

r2/3=i(A, p, pi, P2) = ^i3=i{p, -A, Pi, 



(4.48) 

(4.49) 
(4.50) 
(4.51) 

(4.52) 



due to the p.B ^ Pi symmetry for /3 = 1. Note that this symmetry has nothing to do with 
the orbifold equivalence. 



'For p 7^ 0, the partition functions are not identical even in the leading exponential behavior. 
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4.3.4 Nonperturbative orbifold equivalence between /3 = 4, /3 = 2, and (3 = 1 

Comparing ( 4.391 ), ( 4.42| ), and ( [4.481 ), and using the fls ^ Ui symmetry for /? = 1 RMT, 



one finds tlie relation (note that A = Oat/UB = and p = at /U/ = 0): 

r2/3=4(/iB)U/=0 = 20/3=2 (/^/)Ui3=0 = 2^2/3=1 (/is) |^j=o, (4.53) 
for 2(7/3=4 = (J/3=2 = Cr/3=1, 2A/3=4 = /9/3=2 = A/3=i . (4.54) 

Unhke the relation ( [4 .451 ), this is valid not only for p = (or A = 0) but also for /) 7^ (or 
A 7^ 0), as is consistent with our claim in Sec. |3.3[ The difference of the factor 2 in ( 4.53| ) 



originates from the recipe ( |3.15| ) , or the fact that the /3 = 4 RMT with the size of the matrix 
$ being 2N includes 2 copies of the /3 = 2 or /3 = 1 RMT with the size A^. The factors 2 
in (4.54) originate from the identifications of the chiral condensate and diquark (or pion) 



condensate according to the recipe ( |3.15| ) again, e.g., (nti)/3=4 = 2{uu)p=2 = '2,{uu)p=i. The 
relation ( [4.53| ) leads to the equivalence of the magnitudes of the order parameters (up to 
the factor 2) and the critical temperatures: 

2(7/3=4 (/iB)U/=0 = O-,3=2(A/)Us=0 = 0-/3=l (Ab) U/=0 , (4.55) 
2A/3=4(/iij)|^,=0 = /3/3=2(A/)Us=0 = A/3=i(/iB)|^,=o, (4.56) 

Tp=i{P-B)\fii=o = r^=2(A/)Ufl=o = T'^=i{l^B)\fii=o, (4.57) 

r^=4(AB)U/=0 = 7|'=2(A/)Us=0 = T^=l{P-B)\fli=Q, (4.58) 

as a consequence of the orbifold equivalence. We note that, the equivalence of the neutral 
order parameters and the critical temperatures should be satisfied in the original QCD and 
QCD-like theories as we claimed in Sec. ^ while the effective potentials will not necessarily 



coincide in QCD as (4.53). In the case of the RMT, the effective potential is a function of 
only the neutral order parameters and all the moments are identical; as a result, the effective 
potentials must be identical. In QCD and QCD-like theories, the effective potentials are 
functions of not only the neutral order parameters but also non- neutral ones, so the effective 
potentials will not be identical generally. 

4.4 Brief summary 

We have applied the idea of the orbifold equivalence to the chiral random matrix theories 
(RMTs) and constructed the orbifold projections between the RMTs with different Dyson 
indices /3. The equivalence of the order parameters, both the chiral condensate and di- 
quark (or pion) condensate, has been demonstrated by computing the effective potentials 
of RMTs. 

From the viewpoint of the orbifold projection, the construction of whole class of the 
RMTs at finite ^ can be understood in a unified way (see Fig. ^). However, the (3 = 2 RMT 
at large /i_B is the only exception among these RMTs which cannot be obtained through 
the orbifold projection from the parent RMT because of the spontaneous breaking of the 
projection symmetry. This may be the fundamental reason why an RMT at large which 



reproduces the partition function of the color- flavor locked phase in the e-regime |71| has 
not been constructed yet. 
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5. Conclusion and discussions 



In this paper, we have discussed the universahty of the phase diagrams of QCD and QCD- 
hke theories in the large- A'c hmit via the orbifold equivalence. The whole phase diagrams 
described by the chiral condensate and the superfluid gap are identical between SU(A'c) 
QCD at finite isospin chemical potential and S0(2A'^c) and Sp(2A'^c) gauge theories at 
finite baryon chemical potential fiB- The phase diagrams of these theories outside the BEC- 
BCS crossover regions are also identical to that of SU(A''c) QCD at finite fiB- Especially, 
the chiral condensate and its critical temperature in QCD at finite /Ub should be exactly 
described by those of sign- free QCD at finite fij: the phase-quenched approximation for 
these quantities is exact in the large- A^c limit outside the BEC-BCS crossover region. We 
have also checked that the equivalence is well satisfied for A'^c = 3 at asymptotically high 
densities using the controlled weak-coupling calculations. This leads us to expect that the 
phase-quenched approximation for the chiral condensate also works well even for Nc = 3. 

Our results provide a way to evade the sign problem in the lattice QCD simulation at 
finite fiB, especially in the region relevant to the physics of the chiral phase transition at 
high temperature. The putative QCD critical point may be investigated by studying the 
sign-free QCD at finite /i/ and S0(2A'^c) and Sp(2A'^c) gauge theories at finite hb outside the 
BEC-BCS crossover regions. The lattice QCD simulations in QCD at finite fij were already 
performed in |74, ^] which seem consistent with the results at finite hb though 



they may not be conclusive. Further investigations in this direction would be desirable. 

Other interesting phenomena in QCD at large fj,B and at low T, such as the color 
superconductivity ||^ and the quarkyonic phase |78, 7£], are unfortunately inside the BEC- 



BCS crossover region of the parent SO and Sp gauge theories where the orbifold equivalence 



breaks down. One may add appropriate deformation to S0(2A''c) [^H, ill and Sp(2A''c) 
gauge theories to prevent the BEC of diquark pairing at small fi, which allows us to study 
the properties of QCD beyond fi = m,r/2 using the lattice technique. One should note, 
however, that it is this BEC-BCS crossover region inside which the QCD phase diagram at 
finite hb crucially depends on Nc- For example, the color superconductivity is no longer 
energetically favorable and is replaced by the inhomogeneous chiral density wave in the 
't Hooft hmit [^, ||]. This is in contrast to QCD at finite /is outside the BEC-BCS 
crossover region, QCD at finite fij, and S0(2A''c) and Sp(2A'^c) gauge theories at finite fiB, 
where the phase structures are not affected by Nc qualitatively. There might be some 



^"One might suspect that, for /i < and small but finite T, the system is a gas of baryons in QCD 

at finite fiB, which is completely different from a gas of pions in QCD at finite fii, and the equivalence 
would not hold. However, this difference is irrelevant to the orbifold equivalence of the neutral operators 
in the large- TVc limit. For example, consider the chiral condensate at small T. In QCD at finite /is, the 
thermal excitation of heavy baryons with the mass ~ A^^ is suppressed at small T = 0{Nc) so that the 
chiral condensate remains unchanged from the value in the QCD vacuum. On the other hand, in QCD at 
finite fii, the thermal excitation of pions with the mass ^ cannot change the chiral condensate ~ A'^^. 
Therefore, the magnitude of the chiral condensate should be the same in both theories in this region, and the 
equivalence is satisfied rather trivially. The prediction of the orbifold equivalence is nontrivial at larger T 
near the critical temperature Tc{i-i) of the chiral phase transition; if it is larger than the critical temperature 
of the deconfinement phase transition, it generally depends on /i, for which the orbifold equivalence still 
predicts the exactly the same Tc{fi) in both theories. 
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connection between the region where the physics dramatically changes depending on Nc 
and the region where the orbifold equivalence breaks down. 

As we have revealed in this paper, the idea of the orbifold equivalence is useful to 
discuss the universal properties of different quantum field theories. One should be able to 
see the universality of phase diagrams in QCD and QCD-like theories within a holographic 
model of QCD, the Sakai-Sugimoto model and its generalizations to S0(2A^c) and 
Sp{2Nc) gauge groups |81|. It would also be interesting to argue possible universal proper- 
ties of other systems than QCD. For example, one may generalize the orbifold equivalence 
in the chiral random matrix theories to other class of random matrix theories (Wigner- 
Dyson type and Bogoliubov-de Gennes type) relevant to other systems. 



Note added 

Further study |82] reveals that the QCD critical point is theoretically ruled out in QCD 
at finite hb outside the BEC-BCS crossover regions in the corresponding phase diagrams 
of QCD at finite /i/ and S0(2A^c) and Sp(2A^c) gauge theories at finite /is, at least in the 
large- A'c limit. 
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